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Abstract 

For the frustrated two-dimensional 5 = 1/2 antiferromagnetic Heisenberg 
model close to quantum phase transition we consider the singlet ground states 
retaining both translational and SU(2) symmetry. Besides usually discussed 
checkerboard, spin-liquid and stripe states an unconventional state with two co- 
existing long-range orders appears to be possible at sufficiently large damping of 
spin excitations. The problem is treated in the frames of self-consistent spheri- 
cally symmetric approach. 

PACS74.72.-h, 71.27.+a, 75.20.-g 

1. The general quantum phase transition problem [TJ EJ, [3] is commonly considered 
in the frames of 2D Heisenberg quantum antiferromagnet. The main interest in this 
context is concentrated around the maximally frustrated point J2/J1 ~ 0.5 (frustration 
parameter p = J2HJ1 + J2) ~ 0.3). In this region several energetically competitive 
states with different symmetry appear to be possible ground state ^ atT = 0. That is 
why the recent theoretical efforts are aimed at quantum phase transition near p ~ 0.3 
with several possible scenarios. 

There is firstly a class of states which brake either translational T\ (1 - lattice 
vector) symmetry or spin SU(2) symmetry of the Hamiltonian. These are semiclassical 
Neel checkerboard state (p < 0.3) and semiclassical stripe phase (p > 0.3), where both 
T\ and SU(2) symmetries are broken. Another states of the class are singlet (S*) = 
valence bond crystal (VBC) states (p ~ 0.3, T\ symmetry is broken, SU(2) symmetry 
is restored). Spins in VBC are coupled in pairs forming singlet valence bonds, the 
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latter being arranged in a periodic pattern. Usually only the simplest - columnar and 
plaquette - VBC are considered. The aforenamed states are sketched in Fig.l. 

Let us note, that an accurate comparison of competitive states is of course permis- 
sible only within one and the same method. 

Hereinafter we will consider only another class of pretender states - singlet states 
which differ from semiclassical ones and do not brake neither T\ nor SU(2) symmetry. 
The simplest state of this class can be constructed in the frames of the mentioned 
nearest neighbour valence bond approach. The wavefunction of this state is a sum of 
terms represented in Fig. 2a. Each term is given by disorderly placed valence bonds, 
but the whole spin liquid (SL) state picture is periodical. 

One can generalize the nearest neighbour valence bond concept and consider also 
arbitrary range valence bonds (see Fig.2b). In the latter picture spin-spin correlation 
functions ( ( S^|S , Q) r _ i . 00 for infinitely distant sites can become non-zero under certain 
circumstances, thus opening the way to consider the transition between spin liquid 
singlet state without long-range order (LRO) to singlet state with LRO. 

In the present work we investigate singlet states with and without LRO in the 
frames of Kondo-Yamaji-Shimahara-Takada self-consistent spherically symmetric the- 
ory (SST) [HE], adapted for frustrated Heisenberg model [SI El- SST does not brake 
neither translational nor SU(2) symmetry. The method describes the singlet ground 
state in terms of spin-spin correlation functions and allows to find the triplet exci- 
tations spectrum. Note, that the arising description of Neel and stripe states differs 
from common semiclassical ones, presented in Fig.l, namely, a mean site spin in SST 
is equal to zero (S*) = in any phase, sublattices are absent and the long-range order 
(or its absence) is controlled by spin-spin correlators (S* |S'o) r - > oo for infinitely distant 
sites, or equivalently by the presence or absence of the gap in the spin spectrum at the 
symmetrical points (±7r, ±7r), (0, ±7r), (±7r, 0). 

In more detail, at zero temperature the structure factor c q = (S^St^ in SST is a 
superposition of smooth part c 0q and, for long-range order cases, 5-terms c^(q — q^) 
at symmetric points qQ of the Brilloine zone. This means that in the coordinate space 
the spin-spin correlation function acquires the form 

cr = (S z n+r \S z n ) = 1 J2 co q e iqr + £ c&e**' (l) 

q v 

The first term in the right-hand side of (CQ) defines local correlations and it vanishes 
at infinite distance r — > oo, the second (condensation) part controls LRO. The existence 
of non-zero condensation part is equivalent to zero spin gaps at the appropriate points 
qp of the Brillouin zone. 

In particular for spherically symmetric Neel-type state (SST-Neel) qo =1 = Q = 
(jr, 7r) and the long-range order has a checkerboard pattern 

(S z r \S z ) lrHoo = C d (-ir* +n y (2) 

where r = n x g x + n y g y , and g x = (1, 0), g y = (0, 1) are lattice vectors. The effective 
magnetization m Neel = a/ c§ eel depends on frustration. SST-Neel state is realized at 
low frustration. m Neel decreases from its maximal value at p = to zero at point p* N , 
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where LRO disappears (p* N ~ 0.2, [HIE])- Let us remind, that in SST (S*) = in any 
phase. 

Spherical symmetric stripe state (SST-Stripe) differs from the presented in Fig.l 
semiclassical picture also in the following sense - it is the coherent superposition of 
stripe pictures with horizontal (as in Fig.l) and vertical stripes. For this state qo~ 2,3 = 
X = (0, 7r), (tt, 0) and the long-range order pattern is 

Stripe 

(m)lrHoo = ^-[(-ir + (-1)""] (3) 

SST-Stripe state is realized for higher frustrations. The effective magnetization 
m stnpe _ J c ^ tn P e increases from zero at p* s ~ 0.5 [7J to the maximal value at p = 1. 

The solution with q^ eeZ = Cq*" P6 = 0, i.e. without long-range order, corresponds to 
SL. In this case spin gaps in the symmetrical points of the Brilloine zone are opened, 
their value depending on frustration. The correlation length is defined by the lowest 
spin gap and also depends on frustratrion. In this sense the solution involved is the 
arbitrary range spin liquid (compare Fig. 2b). Self-consistent SL solution is realized in 
the intermediate frustration region p ~ 0.3 E], [9]. 

Note, that besides three aforenamed types of solutions the approach allows to con- 
sider the exotic state with two mutually penetrating long-range orders - checkerboard 
and stripe. The spin gaps in this state must be closed both at Q and at X. The main 
result of this work is that at sufficiently large damping of spin excitations the solution 
of this type does exist and at the point of quantum phase transition is degenerate with 
another (two) solutions. 

The paper is organized as follows. In part 2 we briefly remind the SST calculation 
procedure in the mean field realization. In part 3 the damping of spin excitations is 
introduced in the simplest semiphenomenological way. As it is shown in part 4, at 
sufficiently large but still realistic damping there appears a new type of self-consistent 
solution - with two long-range orders. In passing it is demonstrated that the width of 
spin liquid solution region strongly depends on damping. This part is also devoted to 
brief discussion. 

2. The Hamiltonian of the frustrated Heisenberg model for S = 1/2 spins on a 
square lattice is 

V = E ^ i"S i+g + ~ h J2 S i S i+d (4) 

here the first term in the right-hand side stands for the nearest-neighbour interac- 
tion (Ji - antiferromagnetic exchange constant, g - nearest-neighbour vector), and the 
second - for the next-nearest-neighbour one (with corresponding quantities J2 and d). 
Standard parameter p is a measure of frustration p = J 2 / J, J\ — (1 — p)J, J2 = pJ, 
J = J 1 + J2. All the energetic parameters are expressed in the units of J, hereafter we 
put J = 1. 

In the frames of SST the calculations come to the chain of equations of motion 

00 

for spin retarded Green's function G z nm = (S^\S^) u+ u = -i J dte iujt ([S?(t), S?]). The 
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chain is closed at the second step by the approximation of the following type 



>S'n+gi+g2 , S'n+g 1 'S'n ~ a g($jl ( ^n+gi+g 2 ^'n+g 1 / . . 

/ \ ' ' 



where a r - so called vertex corrections [5] . In the mean- field approximation Green's 
function G z (q, uj) = (S*\St J)^ has the form 



H v q 



where the numerator F q and the spin excitations spectrum u; q are functions of 
lattice sums, involving spin-spin pair correlation functions c r = (S^ +r \S^) for first five 
coordination spheres (see [SI JUl] for details). 

G z (q, uj) defines the structure factor, i.e. Fourier-transform of a correlation function 

c q = (S^SiJ = -- frf W mHlmG 2 (q, W ) = A(2 m ( Wq ) + l) ; ( 7 ) 

m(c<j) = {e% — l) 
which, in turn, leads to five self-consistent equations: 

c r = ^^c q e^ r ; (r = g,d,2g,2d,g + d) (8) 
q 

The system of equations (jSJ) is then solved numerically. 

At T = the structure factor c q ([7]) is a superposition of smooth part and, possibly, 
5-terms Cg<5(q — qo) at symmetric points qo of the Brilloine zone. The latter solution 
corresponds to the presence of long-range order, the condensation part appears in 
Eq.(JS]), and (jSJ) transforms into 

^£ c ° q eiqr + £ c o VqSr (9) 



N 
q 



where co q = defines local correlations, and Cg controls the long-range order. 
Non-zero condensation part leads to additional condition in the self-consistent system 
of equations (JSJ) — the zero spin gap at the appropriate point of the Brillouin zone 
(qp =1 = Q = (-7T, 7r) for SST-Neel and qo~ 2,3 = X = (0, 7r), (it, 0) for SST-Stripe phase). 
Solution with all Cq = 0, i.e. without long-range order corresponds to spin liquid. As 
it is seen from ([9]), in particular, the scheme allows to search for the states with several 
mutually penetrating long-range orders, defined by the points q^. 

3. The described starting approach leads to spin excitations without damping, 
being the mean-field realization of SST. Nevertheless, at this stage it is possible to 
trace the evolution of the ground state with the frustration increase from SST-Neel 
phase to spin liquid and further to SST-Stripe phase [TJ [TUJ . 
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But, as it was shown in [SJ, [TU], the damping of spin excitations is qualitatively 
important, even if it is introduced in the semiphenomenological way. The scheme, 
based on formally exact expression for the Green's function [TT] , in the simplest case 
is reduced to the following Green's function form (instead of (JBJ)): 

G z {^u)= - ^ q (10) 

here only the imaginary part of the polarization operator is taken into account 
and it is written down as the trivial odd function of to. Damping constant 7 must be 
considered as an external parameter, defined by whatever additional arguments. 

The rest of the self-consistent calculations procedure remains the same (of course, 
it must be recalculated at any fixed 7). 

4. As it was mentioned above, in principle one can seek four types of self-consistent 
solutions. The first is the spin-liquid solution without long-range order C| r |^oo — > 0. 
There are also two types of the "sole" long-range order solutions: with checkerboard 
C|r|->oo ~ e iQr ~ (-i)"*+"y and quantum stripe C| r |^oo ~ e iXr ~ ^[(-l) Ux + (-l) ny ] 
correlators pattern. The fourth possible type of solution corresponds to two mutually 
penetrating long-range orders, when the spin-spin correlator at infinitely distant sites 
is the linear combination on checkerboard and stripe laws. 

In the mean-field approximation with zero damping only three self-consistent solu- 
tions are realized in SST - SST-Neel for small frustration p, spin-liquid for intermediate 
p and SST-Stripe for large frustration [7]. 

Fig. 3 represents the energy for all types of self-consistent solutions (found in the 
present work) in the whole range of frustration parameter p for different values of 
damping parameter 7. It is seen, that at small 7 (7 < 0.3) the picture remains 
qualitatively unchanged: for different values of 7 only three types of solutions are 
realizes. 

With growing damping 7 the points of transitions p* N (Neeb — ^Liquid) and p* s 
(Liquids — ^Stripe) are shifted towards each other. It was shown in our previous work 
[8j, that p-region of SL solution shrinks from the left side with the increase of 7, i.e. 
the value p* N increases. Current results show, that the same is true for the right side 
of SL region, corresponding to SL< — ^Stripe transition. 

But at 7 > 0.3 for intermediate frustrations p ~ 0.3 there appears another self- 
consistent solution - with two coexisting long-range orders. For 0.3 < 7 < 0.6 this 
solution is metastable in the whole region of its existence and is disconnected by energy 
with the other solutions. 

At 7 > 0.6 this biordered solution is energetically connected with the SL one for 
p < 0.3 and with stripe state for p > 0.3, in the intermediate interval it remains 
metastable. This is clearly seen from Fig.4, where for 7 = 0.65 the energies of all four 
solutions - SST-Neel, SST-Stripe, SL and biordered are depicted. The corresponding 
condensates for long-range ordered solutions are also shown. 

Note that the value 7 ~ 0.6 can be considered as absolutely realistic [T2], [13]. 

As it is seen from Figs. 3-4, three states of different types (SST-Neel, SL and 
biordered) at 7 > 0.5 are degenerate at the quantum transition point p QPT w 0.28. 
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Thus we come to a conclusion, that in the maximally frustrated region quantum fluc- 
tuations in addition to usually discussed states (checkerboard, SL and other mentioned 
above states) can select the exotic state with two long-range orders. 

Let us mention, that in the biordered state at p > pQ pT with decrease of p the 
checkerboard condensate increases and the stripe condensate rapidly decreases. 

The typical spin excitations spectrum u; q for this state is presented in Fig.5. The 
spectrum is gapless at points X and Q. 

As it was noted above for intermediate damping 0.3 < 7 < 0.6 the biordered 
solution is disconnected by energy with other ones (see Fig. 3, curves for 7 = 0.35 
and 7 = 0.45). To reveal the states which could be energetically connected with 
the biordered state, one must examine other candidates, for example the helicoidal 
solution with qQ = (q, tt), (vr,g). Its semiclassical prototype is the helicoidal state, 
which is realized in Ji — J 2 — J3 Heisenberg model [Hj. Nevertheless this speculation 
is a matter of further study. 

In conclusion, we have found that in the frames of self-consistent spherically sym- 
metric approach for the 2D frustrated Heisenberg model at T = the unconventional 
state with two coexisting long-range orders can appear near the point of quantum 
phase transition p® PT ph 0.28 at sufficiently large damping of spin excitations. 

The work was supported by the Russian Fund of Fundamental Investigations. 
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Figure 1: Commonly considered states for 2D Heisenberg model. Upper left - semi- 
classical Neel checkerboard state, upper right - semiclassical stripe state, thick arrow 
implies the transition with the frustration increase. Lower left - columnar VBC, lower 
right - plaquette VBC 
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Figure 2: a) Nearest neighbour valence bond spin-liquid state, b) Arbitrary range 
valence bond spin-liquid state. The oriented bond between two sites i and j stands for 
singlet valence bond (1/ v / 2)(| hij) - | lit?))- 
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Figure 3: The energy for all found types of self-consistent solutions for < p < 1 and 
different values of damping 7. Solid lines on the left correspond to SST-Neel solutions, 
solid lines on the right - SST-Stripe, intermediate dotted lines are SL solutions, small 
solid curves above SL at p ~ 0.3 correspond to biordered solution. The ordinate axis 
scale fits the curves with 7 = 0.001, other lines are downshifted by O.27. 
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Figure 4: The energy of SST-Neel, SL, biordered and SST-Stripe states in the vicinity 
of quantum phase transition for damping 7 = 0.65. 

The upper graphs show the condensates for the LRO solutions. For checkerboard LRO 
C| rHoo ~ C^(-l) n * +n y, for stripe LRO C| rHoo ~ C s (l/2)[(-l) n - + (-1)""]. Solid lines 
are for the SST-Neel (left) and SST-Stripe (right) states. The dashed and dashed- 
dotted curves are checkerboard and stripe condensates for the biordered state. 



10 



1.2 -, 

co(q)/J 




Figure 5: Spin excitations spectrum a; q for the state with two coexisting long-range or- 
ders. Frustration and damping parameters are p = 0.285 and 7 = 0.65. The spectrum 
is gapless at points X and Q 
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